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We use lattice QCD to fully determine the CKM matrix. \V c d\> \V CS \, \V u b\, \V c t,\ and \V US \ are, 
respectively, directly determined with recent lattice results for form factors of semileptonic 
D — > -nlv, D — > Klv, B —* nil/, B — > Dlv and K — > nlv decays obtained by the Fermilab 
Lattice, MILC, and HPQCD Collaborations. In addition, \V u d\, \Vtb\, \Vts\ and |Vtd| are 
determined by using unitarity of the CKM matrix and the experimental result for sin (2/3). 
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1 Introduction 

The Cabibbo-Kobayashi-Maskawa (CKM) matrix, which relates the mass eigenstates and the 
weak eigenstates in the Standard Model electroweak theory, is a set of parameters. To deter- 
mine each CKM matrix element, one requires both theoretical and experimental inputs. On the 
theoretical side, one needs to know relevant hadronic amplitudes, which often contain nonper- 
turbative QCD effects. A major role of lattice QCD is to calculate such hadronic amplitudes 
reliably and accurately, from first principles. One can then extract the CKM matrix elements 
by combining lattice QCD as the theoretical input with the experimental input such as decay 
rates. In this paper, we show that it is now possible to fully determine the CKM matrix, for the 
first time, using lattice QCD. The result for the full CKM matrix with lattice QCD is: 

/ \v ud \ \v us \ \v ub \ \ 

0.9744(5)(3) 0.225(2)(1) 3.5(5)(5) x KT 3 

\Vcd\ \V CS \ \V d} \ 

0.24(3)(2) 0.97(10)(2) 3.9(f)(3) x 1(T 2 

\Vtd\ \Vu\ \Vtb\ 

\ 8.1(2.7) xl(T 3 3.8(4)(3)xl0~ 2 0.9992(0)(1) J 

where the first errors are from lattice calculations and the second are experimental, except the 
one for \Vtd\ which is a combined lattice and experimental error. The results for the Wolfenstein 
parameters with lattice QCD are: 

A= 0.225(2)(1), A = 0.77(2) (7), p = 0.16(28), rj = 0.36(11). (2) 

To directly determine 5 CKM matrix elements (|V^|, \V CS \, \V u b\, \V c t,\ and |14s|), we use 5 
semileptonic decays (D — > nlu, D — > Klv, B — > irlu, B — > Dlv and K — > irlv), for which the 
techniques for lattice calculations are well established, and thus reliable calculations are possible. 
The accuracy of previous lattice calculations was limited by two large systematic uncertainties 
- the error from the "quenched" approximation (neglect of virtual quark loop effects) and the 
error from the "chiral" extrapolation in light quark mass (mi — ► m u d)- Both led to effect s of 
around 10-20%. Recent work by the Fermilab Lattice, MILC, and HPQCD Collaborations^ 
successfully reduces these two dominant uncertainties. The error from the quenched approxi- 
mation is removed by using the MILC unquenched gauge configurations^, where the effect of 
u, d and s quarks is included (rif = 2 + 1). The error from the chiral extrapolation is greatly 
reduced by using improved staggered quarks. With this improved approach, the accuracy of 
the 5 CKM matrix elements is comparable to that of the Particle Data Group^. The results 
for \V u b\, \V c b\ and \V US \ are preliminary. We then use CKM unitarity to determine the other 4 
CKM matrix elements (IV^I, \Vtb\, \Vt s \ and I n this way, we obtain all 9 CKM matrix 

elements and all the Wolfenstein parameters. The results for the D and B decays have been 
presented in Refs.^4 This work is a part of ongoing project of flavor physi cs with lattice QCD 
by the Fermilab Lattice, MILC, and HPQCD Collaborations; see also" 
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Figure 1: Form factors from BK-based (filled symbols and curves) and non-BK-based (open symbols) analyses 

for D — > nil/ (left panel) and B — > -kIv (right) decays. 

2 5 CKM matrix elements from 5 semileptonic decays 

2.1 D — ► -k(K)Iu, \V cd ^\ and B -> ttIv, \V ub \ 

The differential decay rate dT/dq 2 for the heavy-to- light semileptonic decay H — > P/i/ is pro- 
portional to |^j| 2 |/_|_(g |2 )| 2 , where / + is a form factor of the relevant hadronic amplitude defined 
through 

(P\V»\H) = f + ( q 2 )(p H +p P -Ar + M q 2 )A» 
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Here q = pn — Pp, 
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Ph/^h, P± = Pp — Ev and E = Ep. To 



2 

determine the CKM matrix element |VL-| with the experimental rate J 9 ™* dq 2 (dT/dq 2 ), we 

^min m _ 

calculate / +i o as a function of g 2 . Below we briefly describe the analysis procedure in Refs.''^. 
We first extract the form factors /ii and f±, and interpolate and extrapolate the results for 

/|| and f± to common values of E using the parametrization of Becirevic and Kaidalov (BK)E1. 
We then perform the chiral extrapolation (mi — > m ul j) at each E using the NLO correction in 
staggered chiral perturbation theory^. Finally we convert the results for f± and f\\ at mi = m u d, 
to /+ and /q. The results for /+ and /o are parameterized with the BK formal, 



(1 - q 2 )(l - aq 2 



fo(q 2 ) 



fi 



l-<P/(3> 
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where q 2 = q 2 /m 2 H ,. We obtain DEI = 0.64(3), oP^ = 0.44(4), f3 D ^ n = 1.41(6) 

and f^ K = 0.73(3), a D ^ K = 0.50(4), fi D ^ K = 1.31(7) for the D decays, and = 
0.23(2), a B ^^ = 0.63(5), f3 B ^ n = 1.18(5) for the B decay, where the errors are statistical only. 
To estimate the error from the BK parameterization, i.e., the error for q 2 dependence, we also 
make an alternative analysis, where we perform a 2-dimensional polynomial fit in (m^ E(q 2 )) . 
A comparison between two ana lyses is shown in Fig. ^ The results for D decays agree well with 
recent experimental results We then determine th e CKM matrix elements by integrating 
\f+(q 2 )\ 2 over q 2 and using experimental decay rates |4J1 1 | 12 | p or \y ub \ we use a combined 
average of the decay rate for q 2 > 16 GeV 2 in Refs.^^ an( }LL2l We obtain 



\Vcd\ 

from the D decay, and 



0.239(10)(24)(20) , \V CS \ = 0.969(39) (94) (24) 
\V ub \ x 10 3 = 3.48(29)(38)(47) 



(5) 



(6) 



from the B decays, where the first errors are statistical, the second systematic, and the third 
are experimental errors from the decay rates. The systematic errors are dominated by the finite 
lattice spacing effects, i.e., the lattice discretization effects; see Tabled The results for the CKM 
matrix elements agree with the Particle Data Group averages^ with a comparable accuracy. 
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Figure 2: m;-dependence of A/ for K — > 7rZf decay, together with results from Refs.^^^l 
2.2 B -> Dlv, \V ch \ and K -> irlv, \V US \ 

The differential decay rate of B — ► Dlv is proportional to the square of |V^,|.F(i<;), where ^(w) 
is the relevant form factor and w = v ■ v' with v = pb /ub, v 1 = PD/mn- To ex trac t \V c b\, we 
calculate the form factor at w = 1, .F(l), by employing the double ratio method^l The light 
quark mass dependence of J-(l) is mild, and by extrapolating the result linearly to mi — * 
we obtain^ J :B ^ D (l) = 1.074(18)(16), where the first error is sta tistical, and the second is 
systematic. Combining with the experimental result for | V^,| ^-"(1)^^, we obtain 



\V cb \ x 10 2 = 3.91 (07) (06) (34), 



(7) 



where the first two errors from lattice calculation are smaller than the experimental one (third). 

Finally we study the K — > ttIv decay to determine \V US \. The expression for the K — > it 
decay amplitude is given in an analogous way to Eq. @ . We calculate the K° — > ir~ form factor 
at q 2 = 0, /+(0) = /o(0), by employing the three-steps method, as in Ref.^^. To perform the 
mi -^m u d extrapolation for /+(0), we subtract the leading logarithmic correction fi in chiral 
perturbation theory, i.e., define A/ = 1 + fi — /+(0). We make a fit to A/ adopting an ansatz, 
A/ = (A + Bmi)(m s — mi) 2 , where A, B are fit parameters. The mj-dependence of A f a nd the 
extrapolated result are shown in Fig. [3 together with a recent quenched lattice result'^ and an 
earlier result by Leutwyler and Roos^^l Our preliminary res ult^ is A/ = 0.015(6)(9), giving 
j-K°^TT _ 0.962(6) (9), which agrees well with those of Refs.^^l Combining with a recent 
experimental result for |y us |/ + (0)^l, we obtain 



\V US \ = 0.2250(14)(20)(12). 



(8) 



Table 1: Systematic errors in lattice calculations. For comparison, the error for each CKM matrix element by the 

Particle Data Group - is shown in the last row. 



semileptonic decay 


D -> ir{K)lv 


B -> irlv B- 


-> Dlv 


K — > ttIv 


CKM matrix element 


|Kd(s) 


\v ub \ 


\v*\ 


| Vus | 


q 2 dependence 


2% 


4% 




<1% 


mi — > m U( i extrapolation 


3% (2%) 


4% 


1% 


1% 


operator matching 


<1% 


1% 


1% 


<1% 


discretization effects 


9% 


9% 


<1% 




total systematic error 


10% 


11% 


2% 


1% 


error in PDG average 


5%(1%) 


13% 


4% 


1% 



3 Other 4 CKM matrix elements using unitarity and Wolfenstein parameters 



Having the 5 CKM matrix elements directly determined from the 5 semileptonic decays, we can 
check unitarity of the second row of the CKM matrix. Using Eqs. (JSJ) and (J7J), we get 

{\V cd \ 2 + \V CS \ 2 + iKfel 2 ) 1 / 2 = 1.00(10)(2), (9) 

which is consistent with unitarity. Hereafter the first error is from the lattice calculation and 
the second is experimental, unless otherwise stated. 

We now use unitarity of the CKM matrix to determine the other 4 CKM matrix elements. 



\V u d\, \Vtb\ and \ Vt s \ are easily determined: 

\V ud \ = (l-|Ks| 2 -|Kb| 2 ) 1/2 = 0.9744(5)(3), (10) 

|^t6| = {l-\V ub \ 2 -\V cb \ 2 ) 1/2 = 0.9992(0)(1), (11) 

\Vu\ = \V: s V ub + / \V tb \ ~ \V c * s V cb \ / \V tb \ = 3.8(4)(3)xl0- 2 . (12) 

Eqs. (jSJ), (O and © give some of the Wolfenstein parameters, 

A = \V US \ = 0.225(2)(1), (13) 

A = \V cb \/X 2 = 0.77(2)(7), (14) 

(p 2 + V 2 ) 1 / 2 = \V ub \/(AX 3 ) = 0.40(6)(6). (15) 

To extract \Vtd\ and (p, r/), we use the experimental result for sin(2/3) from B — > (cc)K^ decays. 
Prom a unitary triangle analysis with sin(2/3) = 0.726(37)^1 and Eq. (j!5j) . we obtain 

p = 0.16(28), rj = 0.36(11), (16) 

\Vtd\ = 8.1(2.7) xHT 3 (17) 



with a combined lattice and experimental error, completing the full CKM matrix. 
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